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We  classify  the  equilibrium  solutions  of  the  Smoluchowski  equation  for  dipolar  (ex¬ 
tended)  rigid  nematic  polymers  under  imposed  elongational  flow.  The  Smoluchowski 
equation  couples  the  Maier-Saupe  short-range  interaction,  dipole-dipole  interaction,  and 
an  external  elongational  flow.  We  show  that  all  stable  equilibria  of  rigid,  dipolar  rod  dis¬ 
persions  under  imposed  uniaxial  elongational  flow  field  are  axisymmetric.  This  finding  of 
axisymmetry  significantly  simplifies  any  procedure  of  obtaining  experimentally  observ¬ 
able  equilibria. 

Copyright  ©  2007  Hong  Zhou  et  al.  This  is  an  open  access  article  distributed  under  the 
Creative  Commons  Attribution  License,  which  permits  unrestricted  use,  distribution, 
and  reproduction  in  any  medium,  provided  the  original  work  is  properly  cited. 

1.  Introduction 

Nematic  liquid  crystal  polymers  are  viscoelastic  anisotropic  materials  that  have  many  im¬ 
portant  applications  [  1  ] .  The  dynamic  behavior  of  rigid  rod  nematic  liquid  crystal  poly¬ 
mers  is  modeled  by  the  Smoluchowski  equation  [2-4] .  Analytical  results  on  pure  nematic 
equilibria  have  been  obtained  in  a  series  of  papers  [5-11].  In  [12]  a  2D  Smoluchowski 
equation  under  weak  shear  is  analyzed.  In  [13]  a  coplanar  magnetic  field  is  coupled  in 
the  Smoluchowski  equation  to  investigate  the  monodomain  dynamics  for  rigid  rod  and 
platelet  suspensions.  However,  dipole-dipole  interaction  is  not  included  in  [13].  Recently, 
the  authors  [14, 15]  studied  kinetic  equilibria  of  rigid,  dipolar  rod  ensembles  for  coupled 
dipole-dipole  and  Maier-Saupe  short-range  potentials.  This  work  is  a  natural  extension  of 
[14,  15]  to  include  an  external  elongational  flow.  For  completeness  and  for  reader’s  con¬ 
venience,  we  include  aU  the  lemmas  and  theorems  necessary  for  reaching  the  conclusions 
even  though  some  of  the  lemmas  and  theorems  have  appeared  in  previous  works. 

This  paper  is  organized  as  follows.  We  first  briefly  give  the  mathematical  formula¬ 
tion  of  the  Smoluchowski  equation  for  extended  (polar)  nematics  in  Section  2.  Our  main 
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2  Abstract  and  Applied  Analysis 


theoretical  results  of  the  equilibrium  solutions  of  the  Smoluchowski  equation  are  pre¬ 
sented  in  Sections  3,  4,  and  5.  More  specifically,  we  first  show  that  the  first  moment  of 
an  equilibrium  solution  must  be  aligned  with  one  of  the  principal  axes  of  the  second 
moment.  Then  in  Section  4  we  exploit  free  energy  to  show  that  an  equilibrium  solution 
whose  first  moment  is  not  parallel  to  the  imposed  external  elongational  flow  field  is  un¬ 
stable.  Finally  in  Section  5  we  prove  the  most  important  result  of  this  paper:  all  stable 
equilibrium  solutions  are  axisymmetric.  We  give  concluding  remarks  in  Section  6. 


2.  The  Smoluchowski  equation  for  extended  (polar)  nematics 

In  this  paper,  we  study  equilibrium  solutions  of  the  Smoluchowski  equation  for  rigid 
extended  (polar)  nematics  under  elongational  flow.  Here  the  nematic  molecules  are  mag¬ 
netically  polar.  For  these  extended  nematics,  the  molecular  interaction  includes  both  the 
dipole-dipole  interaction  and  the  short-range  Maier-Saupe  interaction.  We  assume  the 
system  has  an  imposed  elongational  flow  field.  The  potential  due  to  the  external  elonga¬ 
tional  field  is  given  by  [  16] , 


Ve  = -—kTEE:mm,  (2.1) 

where  E  is  the  direction  of  the  elongation,  ao  >  0  corresponds  to  the  stretching  in  the  E  di¬ 
rection  {uniaxial  elongation)  and  cto  <  0  corresponds  to  compressing,  k  is  the  Boltzmann 
constant,  and  T  is  the  absolute  temperature. 

Now  let  p(m,  t)  denote  the  probability  density  function  (pdf)  for  dipolar  rod-like  ne¬ 
matic  molecules  in  unit  direction  m  at  time  t.  The  dynamic  evolution  of  the  orientational 
pdf  for  the  ensembles  of  rigid  rods  with  inherent  dipoles  is  governed  by  the  Smoluchowski 
equation  [3]: 


dt  dm  [kTdm^  dmj 


(2.2) 


Here  d/dm  represents  the  orientational  gradient  operator  [17],  and  the  total  potential  is 
given  by 

17(m)  =  —akT (m)  ■  m  -  bkT (mm)  :  mm  -  ^kTEE  :  mm,  (2.3) 

where  a  is  the  strength  of  the  dipole-dipole  interaction,  b  denotes  the  strength  of  the 
Maier-Saupe  short-range  interaction,  and 

((•))  =  C  {•)pdm  (2.4) 

is  the  ensemble  average  with  respect  to  the  pdf  p,  which  is  a  solution  of  the  Smoluchowski 
equation  (2.2).  For  simplicity,  from  now  on  we  assume  kT  =  1,  or  equivalently  we  assume 
that  all  energies  are  normalized  by  kT. 

The  total  potential  (2.3)  can  be  split  into  the  sum  of  the  internal  potential  representing 
the  interaction  among  polymer  rods  with  themselves  and  the  external  potential: 


17(m)  =  17int(m)  -H  !7ext(m). 


(2.5) 
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where  the  internal  and  external  potentials  are  given  by 

!7ext(m)  =  -  y  EE :  mm, 

!7int(m)  =  !7i(m)  +  L72(m), 

(2.6) 

I7i(m)  =  -a(m)  ■  m, 

1/2(111)  =  -h(mm)  :  mm. 

In  this  paper,  we  restrict  our  study  to  the  case  of  a  >  0,  ao  >  0,  and  h  >  0.  In  particular, 
ao  >0  (uniaxial  elongation)  will  help  us  eliminate  many  unstable  equilibrium  solutions. 


3.  Equilibrium  solutions  for  extended  nematics 

An  equilibrium  solution  of  (2.2)  is  given  by  the  Boltzmann  distribution  [3] 


Peq(iii)  =  -^exp[-  L/(m)],  Z  =  J^exp[  -  !7(m)]dm. 
The  nonlinear  integral  equation  for  the  first  moment  vector  (m)  is 


(m>  =  mpeq(in)dm, 

J  s 

whereas  the  nonlinear  integral  equation  for  the  second  moment  tensor  (mm)  is 


mm)  =  mmpeq(ni)dm. 


(3.1) 


(3.2) 


(3.3) 


We  establish  the  coordinate  system  as  follows.  We  select  the  uniaxial  elongation  field 
E  as  the  z-axis.  We  select  the  x-axis  and  the  y-axis  to  be  perpendicular  to  the  z-axis  but 
otherwise  arbitrary.  In  this  Cartesian  coordinate  system,  we  have 


m  =  (mi, m2, m3),  E  =  (0,0,1),  (m)  =  iqi,q2,q3), 


t/(m)  =  -a{qim\  +  q2m2  +  q^tn^)  -  h(mm)  : 


ao  2 
mm  — —mi 


(3.4) 


The  nonlinear  equations  for  the  first  and  second  moments  become 

qi  =  (mi)  =  miPeq(m)dm, 

J  s 

(mm)i^-  =  (mj-mj)  =  J  mimjpi.q{m)dm. 


(3.5) 


We  select  the  spherical  coordinate  system  {(j),d)  using  the  y-axis  as  the  pole.  In  this 
spherical  coordinate  system,  we  have 


(mi, m2, m3)  =  (sin<jf)sin0,cos0,sin(/)cos0). 


(3.6) 


4  Abstract  and  Applied  Analysis 

Theorem  3.1.  For  an  equilibrium  probability  density,  L/exi  satisfies 

=0.  (3.7) 

Proof.  We  first  prove  that 

^^L7(0,0))  =0.  (3.8) 

This  is  not  surprising.  Physically,  this  quantity  is  the  (negative)  total  torque  on  the  sys¬ 
tem  about  the  y-axis.  Since  the  system  is  in  equilibrium,  the  total  torque  should  be  zero. 
Mathematically,  we  have 


^  ilo  lo  U {(f), d)]dd  sin (l)d(f> 

1  r) 

“v  ^exp[  -  !7((jf),0)]d0sin(/)d(/)  =  0. 

Z  Jo  Jo  ou 


The  second  part  is  to  show  that 


(3.10) 


Again,  this  is  not  surprising  at  all.  Physically,  this  quantity  is  the  torque  on  the  system 
from  the  mutual  interaction.  If  there  is  a  torque  on  polymer  rod  A  from  other  polymer 
rods,  then  by  Newton’s  third  law,  rod  A  exerts  a  torque  of  the  opposite  sign  on  other  rods. 
Thus,  the  sum  of  the  torques  on  all  rods  from  the  mutual  interaction  is  zero.  Mathemati¬ 
cally,  we  prove  it  as  follows.  We  notice  that 


dmi 


(3.11) 


Using  this  fact,  we  immediately  obtain 


=  4  [  -cciqimi  -  q^mi)  exp  [  -  U(m)]dm 
Z  Js 


=  -aiqimi  -  q-^mf)  =  -aiqiq^  -  q^qi)  =  0. 
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Here  we  have  adopted  the  more  concise  and  more  convenient  way  of  writing  the  spherical 
integrals  in  terms  of  m  instead  of  For  the  second  part  of  the  internal  potential,  we 
have 


^  “  ld(m)]dm  =  (mm>  : 


2 

mf 


mim2  mim^X 
=  Itriimi  m\  mama  \  : 
\m1m3  mama  m\  / 


I  2mima 

mama 
\ml  -  m\ 


Combining  (3.8)  and  (3.10)  leads  immediately  to  (3.7). 


m2mj, 

0 

-m\m2 


-mima  )  =  0. 
-Imitni/ 

(3.13) 

□ 


Theorem  3.2.  For  an  equilibrium  probability  density,  the  z-axis  is  an  eigenvector  of  the 
second  moment  tensor  (mm). 


Proof  Differentiating  the  external  potential,  we  have 


de 


Hext(0,d) 


de 


=  aomima. 


(3.14) 


Applying  Theorem  3.1,  we  obtain 


(mima)=0.  (3.15) 

Since  the  x-axis  and  the  y-axis  are  selected  arbitrarily,  exchanging  the  roles  of  the  x-axis 
and  the  y-axis  yields  immediately 


(mama)  =  0. 


Thus,  the  second  moment  tensor  (mm)  has  the  form 


(mm) 


/  (ml)  (mima)  0  \ 
(mima)  (ma)  0 

V  0  0  {ml)l 


(3.16) 


(3.17) 


It  is  obvious  that  the  z-axis  is  an  eigenvector  of  the  the  second  moment  tensor.  □ 

To  facilitate  the  analysis,  let  us  select  the  x-axis  and  the  y-axis  properly  in  the  xy- 
subspace  to  diagonalize  the  matrix  (mm).  If  (m)  is  not  parallel  to  the  z-axis,  then  we  can 
select  the  positive  directions  of  x-axis  and  y-axis  such  that  qi  >0  and  q2  >  0.  In  the  new 
Cartesian  coordinate  system,  we  have 

[/(m)  =  -a{qimi+q2m2  +  q3mi)  -  b(sim\  +  S2ml+  Ssml)  -  ym^  (3.18) 


where 


Sj  =  {mj). 


(3.19) 


6  Abstract  and  Applied  Analysis 


Our  next  move  is  to  prove  two  theorems  which  establish  the  relationship  between  the 
first  moment  vector  and  the  second  moment  tensor. 


Theorem  3.3.  If  an  equilibrium  solution  satisfies  qi  >  0,  then  q^  =  0. 

Proof  We  prove  the  theorem  by  contradiction.  Suppose  there  is  a  solution  satisfying  qi  > 
0  and  q3  f  0.  We  are  going  to  show  that  qi  >0  and  qs  f  0  lead  to  {\/ qf){m\mf)  >  0,  which 
contradicts  with  the  result  of  Theorem  3.2: 


1  If 

—  (mimf)  =  — —  mim3exp[  -  U{mi,m2,mi)]din 
qs  q-iZ  is 

=  -^\  mim3{exp[-U{mi,m2,mi)] 
^3^  J  mi  >0 


exp[  -  U{  -  mi,m2,m2,)]}dm. 

(3.20) 


In  the  above,  the  second  factor  of  the  integrand  is 

^1  (mi, m2, m3)  =  exp[  -  L7(mi,m2,m3)]  -exp[  -  t7(  -  mi, m2, m3)] 

=  2sinh  (aijimi)  exp  (aq^mf)  exp  [aq2m2  -  L/ext(m)  -  L/2(m)]. 

(3.21) 


Notice  that  both  17ext(m)  and  U2(ni)  are  even  functions  of  mi,  m2,  and  m3.  Substituting 
into  (3.20),  we  get 


—  (mim3)  =  ^  f  mim3^(mi,m2,m3)dm 
^3  ^3^  Jml>0 


=  I  mim3{^(mi,m2,m3)  -g{-  m]_,m2,mi)}dm 

sinh  {aq^mf) 


qsZ  .  ^ml>0,m3  >0 

4 

2jml>0,m3>0  ^  qs 

X  exp  [aq2m2  -  f/ext(m)  -  L72(m)]dm. 


mi  m3  sinh  {aq\mi)- 


Since  the  integrand  satisfies 

sinh  (aq^mf) 


mi  m3  sinh  (aqimi)- 


?3 


exp  [aq2m2  -  t/ext(m)  -  L72(m)]  >  0 

for  mi  >  0,  m3  >  0,  qi  >  0,  q3  f  0 


(3.22) 


(3.23) 


we  have  {\/ qfilmxmf)  >  0,  which  contradicts  with  the  result  of  Theorem  3.2.  □ 

Remark  3.4.  Theorem  3.3  tells  us  that  if  the  first  moment  vector  (m)  is  not  parallel  to  the 
z-axis,  then  it  must  be  perpendicular  to  the  z-axis. 

Theorem  3.5.  If  an  equilibrium  solution  satisfies  qi  >  0,  then  q2  =  0. 

Proof  The  proof  of  this  theorem  is  very  similar  to  the  proof  of  Theorem  3.3.  In  the  proof 
of  Theorem  3.3,  we  did  not  use  that  the  z-axis  is  the  direction  of  E.  We  only  used  the  fact 
that  both  [/ext(ni)  and  1/2(111)  are  even  functions  of  mi,  m2,  and  m3.  So  by  exchanging  the 
roles  of  m2  and  m3,  we  can  extend  the  proof  of  Theorem  3.3  to  Theorem  3.5.  □ 
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Remark  3.6.  Theorem  3.5  indicates  that  if  the  first  moment  (m)  is  not  parallel  to  the 
z-axis,  then  it  must  be  parallel  to  either  the  x-axis  or  y-axis. 

Combining  the  results  of  Theorems  3.3  and  3.5,  we  conclude  that  the  first  moment 
(m)  of  the  equilibrium  solution  of  (2.2)  must  be  aligned  with  one  of  the  major  axes  of 
the  second  moment  tensor  (mm). 


4.  Free  energy  and  stability 

In  this  section  we  consider  the  stability  of  an  equilibrium  solution.  The  study  of  stability 
allows  us  to  narrow  down  the  possible  solutions  that  can  be  observed  experimentally.  To 
investigate  the  stability  of  an  equilibrium  solution,  it  is  useful  to  exploit  the  free  energy  of 
the  system. 

To  do  so,  consider  an  arbitrary  probability  density  p(m),  which  is  not  necessarily  an 
equilibrium  probability  density.  The  free  energy  of  the  probability  density  p(m)  can  be 
written  as  [8] 

Glp]  =  J  p(m)lnp(m)dm  -  ^  JJ  m'  ■  mp(m')p(m)dm'dm 

-  -  ff  m'm' :  mmp(m')p(m)dm'dm  -  —  ff  (E  ■  m)^p(m)dm' dm 
2  JJs  2  JJs  (4,1) 

=  f  p(m)lnp(m)dm-  ^(m)  ■  (m>  -  ^(mm)  :  (mm)  -  ^((E  ■  m)^) 

Js  z  z  z 

—  Gent  [p]  +Gi[p]  +  G2[p]  +  Gext  [pl 


In  the  above  expression  Gent[f  ]  corresponds  to  the  entropic  part  of  the  free  energy,  Gi  [p] 
and  G2  Ip]  are  free  energy  parts  associated  with  the  two  mutual  interactions,  and  Gext  [p]  is 
the  free  energy  part  associated  with  the  external  field.  Here  ( ■ )  represents  the  mean  taken 
with  respect  to  the  probability  density  p(m).  For  the  clarity  of  analysis  below,  we  intro¬ 
duce  two  different  notations  for  means  taken  with  respect  to  two  different  probability 
densities. 

(i)  ('leq  denotes  the  mean  taken  with  respect  to  the  equilirium  probability  den¬ 
sity,  whereas  ( ■ )  represents  the  mean  taken  with  respect  to  a  general  probability 
density  (usually  a  perturbed  probability  density  near  the  equilirium  probability 
density). 

Recall  from  Theorem  3.5  that  for  an  equilibrium  solution  of  the  Smoluchowski  equa¬ 
tion  (2.2),  if  the  first  moment  (m)  is  not  parallel  to  the  z-axis  (which  is  the  direction  of 
the  external  elongational  flow),  then  it  must  be  parallel  to  either  the  x-axis  or  y-axis.  We 
call  these  solutions  with  (m)  not  parallel  to  the  z-axis  nonparallel  solutions.  Now  we  show 
that  nonparallel  equilibrium  solutions  are  actually  unstable.  Therefore,  we  can  exclude 
them  from  our  study  and  focus  on  the  parallel  solutions. 

Theorem  4.1.  If  an  equilibrium  solution  satisfies  S3  <  Si  or  S3  <  S2,  then  it  is  unstable. 


Proof  We  present  the  proof  for  the  case  of  S3  <  si.  The  proof  for  the  case  of  S3  <  S2  is 
similar.  □ 
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In  order  to  prove  that  an  equilibrium  solution  with  S3  <  Si  is  unstable,  we  show  that  the 
free  energy  G[p]  does  not  reach  a  local  minimum  at  peq(m).  In  particular,  we  construct  a 
perturbed  probability  density  p(m)  arbitrarily  close  to  the  equilibrium  probability  den¬ 
sity  peq(ni)  such  that 

G[p(m)]  <  G[peq(ni)].  (4.2) 


We  rotate  the  equilibrium  probability  density  peq{m)  about  the  y-axis  by  £  (a  small 
angle)  and  use  the  result  as  the  perturbed  probability  density  p(in).  Mathematically  this  is 
equivalent  to  keeping  the  probability  density  fixed  but  rotating  the  external  elongational 
flow  field  about  the  y-axis  by  —e.  After  the  rotation,  the  external  elongational  flow  field  is 

E  =  (sin  £,0,  cose).  (4.3) 


We  have 

Gent[Peq]j 
=  Gl  [/^eq]j 

=  G2  [peq] , 

=  -^((E-m)2)  =  -^((E-m)^),^ 

2  2  4 

=  -^((misin£  +  m3C0S£)^)j.q 

=  -y(('«i)eq  +  £('«l'«3)eq  +  £^[('«?)eq-  ('«3)eq]  +  -  ■  ■) 

=  Gext[feq]  “  (^1  “  +  '  '  '  ■ 

If  S3  <  Si,  then  for  £  sufficiently  small  we  have  Gext[p]  <  Gext[peq]-  It  follows  that  G[p]  < 
G[peq]  and  the  equilibrium  solution  peq  is  unstable. 

Theorem  4.2.  If  the  first  moment  of  an  equilibrium  solution  satisfies  qi  >  0,  then  the  dipole- 
dipole  interaction  strength  a  is  related  to  the  second  moment  through  the  inequality  asi  >  1. 

Proof  Suppose,  on  the  contrary,  that  asi  <  1. 

The  parameter  qi  has  a  fixed  value  in  the  equilibrium  solution  peq(ni).  Here  we  rename 
it  V  and  treat  it  as  an  independent  variable.  We  consider  the  probability  density 

p(m,v)  =  -^expra(vmi  -l-q2m2-l-q3m3)  -l-h(sim^-l-S2m^-tS3m3)  + 

f  ® 

Z  =  J  exp  |^a(vmi  +  q2m2  +  qT,mf)  +  b{s\m\  +S2m\  +S3ml)  +  y 

For  clarity,  let  ( ■  )v  denote  the  mean  taken  with  respect  to  the  probability  density  p(m,  v). 
Note  that  p{m,v)\y=q^  =  peq(m).  Consider  the  function 


Gen.[T] 

Gl[p] 

G2[p] 

Gext[T] 


F(v)  =  V-  (mi)^. 


(4.6) 
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which  satisfies  F(0)  =  0  and  F{qi)  =  0.  We  are  going  to  show  that  asi  <  1  implies  that 
F'{v)  >  0  for  V  e  (0,qi),  which  contradicts  F(0)  =  F(qi)  =  0.  Hence  asi  >  1.  We  do  it  in 
several  steps. 

Step  1.  Differentiating  with  respect  to  v  yields 

^p(m,v)  =  a{mi  -  (mi)  Jp(m, v), 

^(mi)^  =  a{mi{mi  -  (mi)]), 

F'(v)  =  1  -  a{ml)^  +  a{mi)l. 

To  show  F'(v)  >  0  for  v  G  {0,qi),  we  only  need  to  show  a(ml}v  <  1  for  v  G  {0,qi).  Con¬ 
sider  function  g{v)  =  {m\)y.  Because  of  the  supposition  asi  <  1,  we  have  ccgiqi)  <  1. 
Thus,  to  show  ag{v)  <  1  for  v  G  (0,qi),  we  only  need  to  show^'(v)  >  0  for  v  G  (0,(^i). 
Step  2.  g{v)  satisfies  the  property  that  g'  (vo)  =  0  implies  g"  (vo)  >  0.  To  see  this,  differen¬ 
tiating  ^(v)  with  respect  to  v  yields 

/(v)  =  (mi(mi  -  (mi)J)^.  (4.8) 

Differentiating  a  second  time  respect  to  v,  we  get 

^"(v)  =  {mjimi  -  (m?)^(mi(mi  -  (mi)J)^ 

=  {tnt)^-2{ml)^{mi)^+{ml)^{mi)l-  {ml)l+ {mj) 

=  var(mf)  -2/(v)(mi)^. 

Therefore,  if^'(vo)  =  0,  then  we  have^"(vo)  =  var(mi)  >  0. 

Step  3.  When  v  =  0,  the  probability  density  p(m,0)  is  symmetric  with  respect  to  mi.  So 
we  have 


g'(0)  =  (m^)  I  ^^0  -  (mf)  I  I  v=o  =  0-  (4-10) 

Using  the  result  of  Step  2,  we  conclude  that^'(v)  >  0  for  v  G  (0,qi).  This  leads  immedi¬ 
ately  to  F'(v)  >  0  for  V  G  (0,qi),  which  contradicts  F(0)  =  F(qi)  =  0.  □ 

Theorem  4.3.  If  an  equilibrium  solution  satisfies  qs  =  0  and  asa  >  1,  then  it  is  unstable. 

Proof  This  theorem  does  not  specify  any  condition  on  qi.  Recall  the  method  we  used  in 
selecting  the  axes:  if  the  equilibrium  solution  is  not  parallel  to  E,  then  we  select  the  axes 
to  make  qi  >  0  and  q2  >  0.  Theorem  3.5  tells  us  that  if  qi  >  0,  then  q2  =  0  and  q^  =  0.  So 
we  always  have  q2  =  0.  □ 

We  only  need  to  show  that  the  free  energy  G[p]  does  not  attain  a  local  minimum  at 
Peq(ni).  More  precisely,  we  show  that  there  exists  a  perturbed  probability  density  p(m) 
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arbitrarily  close  to  the  equilibrium  probability  density  peq(ni)  such  that 

G[p(m)]  <  G[peq(ni)].  (4.11) 


We  consider 


p(m)  =  (1  +£m3)peq(m). 


(4.12) 


Since  (waleq  =  qs  =  0,  we  have  jgp{m)dm  =  1,  which  means  p(m)  is  a  probability  density. 
We  calculate  the  four  parts  of  the  free  energy  of  the  perturbed  probability  density  p(m). 

Because  q2  =  0  and  =  0,  the  equilibrium  probability  density  peq(m)  is  symmetric 
with  respect  to  m2  and  m3.  Using  the  Taylor  expansion 

1  T 

(fl  +  Ax)ln(a  +  Ax)  =  alna  +  (Ina  +  l)AxH - (Ax)  +  ■  ■  ■  (4.13) 

2a 


we  have 


Gent  [p  ]  ~  Gent  [Peq]  ^  2  ^^^^eq  ’ 

Gl[p]=Gl[peq]-e^“(m^>eq. 


The  symmetry  of  Peq(m)  gives  us 


(mmm3),q 


/  0  0 
0  0 
\(mimi)^q  0 


0 

0 


Substituting  into  G2[p]  yields 


Clip]  =  G2[pe<i]  - 
Gext  [p  ]  “  Gext  [f^eq]  ■ 


(4.14) 


(4.15) 


(4.16) 


Combining  (4.14)  and  (4.16),  we  obtain 
G[p]  -G[peq] 

^  (4.17) 

5  -£^-53(0:53  -  1). 

If  as3  >  1,  then  for  £  sufficiently  small  we  have  G[p]  <  G[peq]j  which  means  the  equilib¬ 
rium  solution  peq  is  unstable. 


Theorem  4.4.  If  an  equilibrium  solution  satisfies  qi  >  0,  then  it  is  unstable. 
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Proof.  Theorem  3.5  tells  us  that  qi  >0  implies  q2  =  0  and  q^  =  0.  We  discuss  two  cases. 

(i)  Case  1:  S3  <  Si.  Theorem  4.1  tells  us  that  the  equilibrium  solution  is  unstable. 

(ii)  Case  2:  S3  >  Si.  Using  Theorem  4.2,  we  have  asi  >1.  It  follows  that  as^  >1.  Using 
Theorem  4.3,  we  conclude  that  the  equilibrium  solution  is  unstable. 

Therefore,  all  nonparallel  solutions  are  unstable.  In  this  way,  we  have  excluded  all  non¬ 
parallel  solutions  from  our  study.  □ 

5.  All  stable  equilibria  are  axisymmetric 

In  the  previous  section  we  have  concluded  that  all  nonparallel  solutions  are  unstable.  So 
from  now  on  we  only  consider  parallel  solutions.  It  is  well  known  that  for  pure  (nondipo- 
lar)  nematic  rod  ensembles  where  (m)  =  0  [8,  11],  the  stable  equilibria  satisfy  either 
5i  =  S2  =  S3  (isotropic  phase)  or  53  >  Si  =  S2  (prolate  uniaxial  phase)  in  the  selected  coor¬ 
dinate  system.  In  [  15]  we  showed  that  this  is  also  the  case  for  extended  (dipolar)  nematics. 
In  this  section  we  are  going  to  show  that  in  the  presence  of  imposed  uniaxial  elongational 
flow  the  system  still  retains  this  axisymmetry. 

Below  for  the  case  of  extended  nematic  equilibria  in  which  (m)  may  be  nonzero, 
we  will  show  that  if  (m)  f  0  (i.e.,  qi  =  qi  =  0  and  ^3  >  0  by  the  selection  of  the  co¬ 
ordinate  system  and  by  the  result  of  Theorem  3.1),  then  a  stable  equilibrium  solution 
must  be  uniaxial.  Furthermore,  the  axis  of  symmetry  must  be  the  major  director  (i.e., 
the  eigenvector  of  the  second  moment  corresponding  to  the  largest  eigenvalue).  That  is, 
(ml)  =  (ml)  <  (ml). 

It  should  be  point  out  that  in  ]15]  we  proved  five  lemmas  that  paved  the  roads  to  reach 
axisymmetry.  In  order  to  prove  axisymmetry  for  the  extended  nematics  in  the  presence 
of  an  external  elongational  flow,  we  will  take  full  advantage  of  the  well-established  results 
in  [15]. 

First,  we  recall  that  in  [15]  there  is  no  external  flow  field  and  the  probability  density  is 
given  by 


(5.1) 


(5.2) 

(5.3) 


Under  the  constraints  5i  <  S3  and  S2  <  S3  we  proved  in  [15]  that  there  is  no  equilibrium 
solution  such  that  Si  <  S2.  This  holds  for  all  a  >  0  and  b  >  0.  Since  the  nonexistence  of 
Si  <  $2  (where  Si  <  S3  and  S2  <  S3)  is  true  for  all  a  >  0  and  h  >  0,  we  introduce  A  =  aq^  as  a 
parameter  and  treat  a  as  unknown.  Equation  (5.3)  becomes 


(5.4) 
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So  (5.3)  can  be  satisfied  by  selecting  a  suitable  value  of  a.  Therefore,  (5.2)  cannot  be 
satisfied  with  5i  <  S2  (where  Si  <  S3  and  S2  <  S3).  Again,  this  is  true  for  all  A  and  b  >  0. 

Let  us  introduce  rj  =  hsj  as  unknowns.  Notice  that  the  pdf  in  (5.1)  does  not  depend 
on  b  any  more  once  we  know  (ri,r2,r3).  Equation  (5.2)  yields 

{m\)  =  {ml)  =  {ml)  =  (5.5) 

It  follows  that  (5.5)  cannot  he  satisfied  with  ri  <  r2  (where  ri  <  r3  and  ri  <  r^).  This  is  true 
for  all  A  and  b  >0.  Further,  we  introduce 

rji=ri-ri<  0,  t]2  =  r2-r^<  0.  (5.6) 

Using  the  fact  that  ml  +  ml  +  ml  =  1,  the  pdf  in  (5.1)  can  be  rewritten  as 

Pi(m)  =  ^exp{Xm3  +  t]iml  +  t]2ml), 

r  (5.7) 

Z=J  exp[Xm3  +  Tjiml  +  rj2ml]dm. 

Note  that  the  parameter  b  does  not  appear  in  (5.7).  Similarly,  (5.5)  turns  into 

{m\-ml)  =  ^,  {ml-ml)  =  ^.  (5.8) 

One  concludes  that  (5.8)  cannot  be  satisfied  with  rji  <  rj2  <  0.  This  conclusion  holds  for 
all  A  and  b  >  0.  Next,  we  prove  two  theorems  for  the  extended  nematics  without  external 
fields. 

Theorem  5.1.  In  the  region  rji  <rj2  <  0,  it  is  true  for  all  A  that 

rj2{m\- ml)  -  tji{ml- ml)  f  Q.  (5.9) 


Proof  Suppose,  on  the  contrary,  that  rj2{ml  -  ml)  -  r\i{ml 

{ml- ml)  {ml- ml) 


ml)  =0.  Then  we  get 

(5.10) 


’ll  m 

Select  the  value  in  (5.10)  as  l/b  and  we  obtain  (5.8),  which  violates  the  fact  that  (5.8) 
cannot  be  satisfied  with  r\\  <r\2<  0. 

Note  that  b  does  not  appear  in  Theorem  5.1.  □ 

Theorem  5.2.  In  the  region  <  f/2  <  0,  it  is  true  for  all  A  that 

t]2{ml  -  ml)  -  rii{ml- ml)  <0.  (5.11) 


Proof  Let  us  denote  t]2{ml  -  ml)  -  rjiiml  -  ml)  by  H{rji,tj2).  It  is  easy  to  verify  that 
fd  ( >7 1  >  >72 )  is  a  continuous  function  of  ( >7 1 ,  >72 ) .  From  Theorem  5 . 1 ,  H  ( >7 1 ,  >72 )  is  nonzero  in 
the  region  >71  <  >72  <  0  and  thus  it  does  not  change  its  sign  in  this  region. 


Hong  Zhou  et  al.  13 

Consider  the  case  where  rji  <  rj 2  and  both  tji  and  tj2  approach  -00.  Then  (m3)  ^  1, 
(ml)  ^0,  and  (ml)  ^0.  Therefore, 


H{rji,ri2)  =  (rji  -  rj2){ml)  +ri2{ml)  -  rjiiml)  <  0  (5.12) 

as  rji  ^  -00  and  >/2  ^  (^1  <  '72)-  This  completes  the  proof  of  Theorem  5.2.  □ 

Now  we  use  the  above  results  to  consider  the  case  of  ao  >  0,  which  corresponds  to  the 
coupling  of  an  imposed  uniaxial  elongational  flow.  In  this  case,  the  pdf  is  given  by 


P2(m)  =  ^exp|^a(j3m3  +  h(simf +52m^+S3m|)  + 

Z  =  J  exp  j^aq3m3  +  h{s\m\  +S2m\  +S3ml)  +  y  m^jdm. 
As  before,  we  introduce 

,  .  OCq  OCq 

A  =  aq3,  rj=bsj,  f]\=ri-r3-—,  >72  =  ?'2  -  r3  -  y . 

Then  the  pdf  (5.13)  becomes 


(5.13) 


(5.14) 


P2(m)  =  -exp[Am3  +  >7imi +?72m2], 

Z=  exp  [Am3  +  f/imj  +  mm^ldm. 

Js 


(5.15) 


Note  that  this  pdf  is  exactly  the  same  as  the  previous  case  (5.7)  with  qi,  ^2  defined  slightly 
differently.  Consequently,  (5.8)  is  also  modified  slightly. 


{ml -ml) 


+  (Xq/2 

b 


{ml- ml) 


r]2  +  oJo/2 
b 


(5.16) 


Our  next  theorem  generalizes  earlier  conclusions  to  include  the  coupling  of  external 
elongational  flow. 

Theorem  5.3.  In  the  region  qi  <  q2  <  0,  it  is  true  for  all  A  and  b  >0  that  (5.16)  cannot  be 
satisfied. 

Proof  We  use  proof  by  contradiction.  Suppose  (5.16)  can  be  satisfied  with  ?7i  <  >72  <  0  for 
some  value  of  b.  Then 


(^2  +  y  )  {ml  -  ml)  -  (^qi  +  y  j  {ml  -  ml)  =  0,  (5.17) 

or 

q2{ml  -  ml)  -  qi{ml  -  ml)  +  ^{ml  -  ml)  =  0.  (5.18) 


From  Theorem  5.2,  q2{ml  -  ml)  -  qiiml  -  ml)  <  0.  Therefore,  {ml  -  ml)  >  0. 
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Next,  we  show  {ml  -  ml)  <0  which  contradicts  the  above  result.  To  do  so,  we  rewrite 
the  pdf  (5.15)  as 


P2(ni)  =  ^exp|^. 


Ama  +  {ml  +  mj)  +  ^  {ml  -  m^)  j , 


and  then 


{ml-m\)  =  - 


=  ^  J  exp  |^Am3+  {ml+ml)~j  {ml -ml)  expj^'^^  {ml-ml) 


(5.19) 


dm. 

(5.20) 


Switching  the  role  of  integration  variables  m\  and  m2  in  the  integral  on  the  right-hand 
side  only  (note  that  mi,  m2  on  the  left-hand  side  are  random  variables  and  have  different 
meanings)  yields 


I  I^Ama-t  {ml+ml)J  {ml-ml)  exp  {ml-ml) 


(ml-ml)  =  2 

Adding  (5.20)  and  (5.21)  and  averaging,  we  obtain 
(ml-ml) 

{ 


dm. 

(5.21) 


= -A  I  exp  I  Am3  +  ^  {ml  +  ml)  (mf  -  m^)  sinh  (wj  -  m^) 


dm  <  0, 
(5.22) 


which  contradicts  earlier  result. 


□ 


6.  Conclusions 

The  stable  equilibrium  solutions  of  rigid,  dipolar  rod  ensembles  (extended  nematics) 
under  imposed  elongational  field  are  shown  to  be  axisymmetric.  Moreover,  the  distin¬ 
guished  axis  of  symmetry  of  stable  anisotropic  equilibria  coincides  with  the  first  moment 
of  the  probability  density  function  (pdf),  the  major  director  of  the  second  moment  of  the 
pdf  (eigenvector  associated  with  the  largest  eigenvalue),  and  the  imposed  elongational 
flow  field.  This  finding  of  axisymmetry  provides  reduction  in  the  degree  of  freedom  in 
the  representation  of  the  pdf  solution  and  thereby  significantly  simplifies  any  process  of 
obtaining  physically  observable  equilibria. 
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